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ABSTRACT 

For studies of Galactic evolution, the accurate characterization of stars in terms of their evolu- 
tionary stage and population membership is of fundamental importance. A standard approach 
relies on extracting this information from stellar evolution models but requires the effective 
temperature, surface gravity, and metallicity of a star obtained by independent means. In pre- 
vious work, we determined accurate effective temperatures and non-LTE \ogg and [Fe/H] 
(NLTE-Opt) for a large sample of metal-poor stars, -3 < [Fe/H] < -0.5, selected from the 
RAVE survey. As a continuation of that work, we derive here their masses, ages, and dis- 
tances using a Bayesian scheme and GARSTEC stellar tracks. For comparison, we also use 
stellar parameters determined from the widely-used ID LTE excitation-ionization balance of 
Fe (LTE-Fe). We find that the latter leads to systematically underestimated stellar ages, by 
10-30%, but overestimated masses and distances. Metal-poor giants suffer from the largest 
fractional distance biases of 70%. Furthermore, we compa re our results with those released 
by th e RAVE collaboration for the stars in common (DR3, IZ witter et al.ll20l61 ISiebert et al.l 
1201 lb . This reveals -400 to +400 K offsets in effective temperature, -0.5 to 1 dex offsets 
in surface gravity, and 10 to 70% in distances. The systematic trends strongly resemble the 
correlation we find between the NLTE-Opt and LTE-Fe parameters, indicating that the RAVE 
DR3 data may be affected by the physical limitations of the ID LTE synthetic spectra. Our 
results bear on any study, where spectrophotometric distances underlie stellar kinematics. In 
particular, they shed new light on the debated controversy about the Galactic halo origin raised 
by the SDSS/SEGUE observations. 

Key words: stars: fundamental parameters — stars: distances — stars: evolution — stars: 
late-type — Galaxy: kinematics and dynamics 



1 INTRODUCTION 

One of the major problems in modern astrophysics is to determine 
physical parameters of stars, masses and ages, and their kinemat- 
ics. Whereas various methods exist, most of them are not suit- 
able for large stellar samples. Thus the most common approach 
is fitting to stellar evolution models, i.e., comparison of an ob- 
served location of a star in the T tS — My plane with isochrones 
or evolutionary tracks for a given metallicity and ^-enhancement. 
This method demands stellar atmospheric p arameters as an input , 
which can be derived fr om photometry (e. g..lHolmberg et al.l2007l ; 
ICasagrande et alfeoill) or spectroscopy dValenti & Fischer||2005I ). 
or so me combination thereof (e.g. application to SDSS/SEGUE 
data in lLee et al]|2008h . The drawback of photometry is that it de- 
pends on reddening, which is difficult to estimate, and demands 
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calibration of metallicity and alpha-enhancement on spectroscopic 
data. Additionally, information on gravity is not obtained directly, 
a quantity needed to accurately determine distances. The alterna- 
tive and preferable method is to use spectroscopy directly, that al- 
lows a consistent determination of all atmospheric quantitites, ef- 
fective temperature, surface gravity, and element abundances, si- 
multaneously. Thus, it minimizes the the total error caused by var- 
ious sources of uncertainties inherent to different methods of pa- 
rameter determination. 

Wit h the advent of larg e-scale stellar surveys such a s 



RAVE l lSteinmetz etail 120061). G CS l lNordstrom etail 120041). 
SDSS /SEGUE tetnnv et al.l |2009|> . GAIA-ESO dGilmore et all 
I2012T) much attention has been devoted to developing sophisticated 
statistical methods to best exploit the available spectroscopic and 
photome tric information. The basic technique of iso chrone fitting 
(see e.g. iNg & Bertellill 19981 : iLiu & ChaboveifeOOOh or a combi- 
nation of complementary techniques jLachaume et al ] fl99i . have 



2 A. Serenelli 



evolved into more systematic methods that can be applied to very 
large amounts of data. Different alt ernatives based on a frequentist 
approach (e.g. lBreddels et al.l2010h or Bayesian analysis of ever in- 
creas i ng sophistication dJergensen & Lindegrenl20 05; Take da et al.l 
120071 : ICasagrande et al.ll201 ll ; iBurnett & Binney||2010l) have been 
resented and their advantages and disadvan tages discussed 
Pont & Evej|2004 |j0rgensen & Lindegrenll2005l) . 

As a result, spectrophotometric determination of masses, ages, 
and distances for late-type stars is nowadays routinely applied in 
large-scale stellar surveys and, although with different implemen- 
tations, lies at the core of almost all studies aimed at characterizing 
stellar populations in the Milky Way, both in the intrinsic stellar 
properties such as age-m etallicity relations dHolmberg et alj|2007l ; 
ICasagrande et al. 1201 1 ) as well as in their kinematical proper - 



ties dCarollo et all20ldTschonrich et al.l201 ltlBurnett et . alj201 ll) . 
However, in these studies little or no attention has been given to the 
physical accuracy of stellar atmospheric parameters, which, in one 
way or another, always rely on models describing radiative trans- 
port in stellar atmospheres. In spectroscopy, all theoretical grids 
of synthetic stellar spectra for late-type stars have been computed 
based on the assumptions of local thermodynamic and ID hydro- 
static equilibrium (hereafter ID LTE). These are, fo r exam ple, the 
widely-used synthetic spectra grids bv lMunari et al.ld2005l) . which 
rely on the Kurucz ID LTE model atmosphere and radiative trans- 
fer codes. However, recent observational and theoretial studies re- 
vealed substantial systematic biases in basic stellar parameters at 
low metallicity and/or gravity caused by the breakd own of LTE and 
ID hydrostatic equilib rium approximations (e.g., lAsplund|[2003 ; 
iBergemann et al.l2012h . It must be pointed out that these biases are 
ubiquitous, and affect the results obtained from the analysis of high- 
and low-resolution spectra. The reason is that departures from ID 
HE and LTE often change equivalent widths (or, equivalently, line 
indices), and consequently abundances inferred from a diagnostic 
spectral feature. This effect does not vanish with a decreasing spec- 
tral resolution, which merely alters the line profile shape without 
affecting the total energy absorbed in a line. Although no detailed 
quantitative estimates have been performed so far, one may also ex- 
pect a non-negligible impact on certain photometric indices. As a 
consequence, it is completely unknown whether, and how, system- 
atic errors in basic stellar parameters caused by the assumptions of 
ID LTE impact the inferred stellar properties, that is, in our context, 
ages, masses, and distance s. 

In lRuchtietalJ l [20ll , hereafter Paper 1), we presented a com- 
prehensive study of systematic errors in the determination of spec- 
troscopic stellar parameters arising from using the classical method 
of excitation-ionization balance of Fe I and Fe II lines based on ID 
LTE model atmospheres. Based on the results for the large sample 
of metal-poor stars (-3 < [Fe/H] < -0.5) selected from the RAVE 
catalogue, we concluded that the impact of using inaccurate effec- 
tive temperature scale and neglecting the effects of non-local ther- 
modynamic equilibrium (NLTE) on the determination of gravities 
and metallicities leads to systematic biases of up to 400 K in effec- 
tive temperature, 1.5 dex in surface gravity and 0.5 dex in metallic- 
itjQ In this work, we use results from Paper 1 to analyze the impact 

1 While the methods developed in the Paper 1 are best suited for applica- 
tion to high-resolution spectra, they can be, in principle, extended to lower- 
resolution data if the wavelength coverage is sufficient to measure the H 
lines or accurate T c g can be obtained by independent means, i.e., from 
high-quality photometry with a careful correction for reddening. A more 
reasonable way is to resort to standard methods relying on global \ 2 mini- 
mization, as e.g. in lZwitter et al.| |2008j), using, however,the synthetic grids 



of using physically accurate spectroscopic stellar parameters on the 
determination of their masses, ages, and distances. 

Our paper is structured as follows: Section |2]gives details on 
the observational data and the different techniques used in the de- 
termination of basic spectroscopic parameters. Section|3]is a brief 
description of the grid of stellar evolution models computed for 
this work. In Section|4]we describe the novel Bayesian method we 
apply to obtain masses, ages and, distances, and discuss the lim- 
itations, particularly imposed by the degeneracy of stellar models 
on the observational plane. We present the results in Section[5]and 
compare them with those published by the RAVE collaboration. 
Finally, in Section [6] we summarize the impact of our findings in 
the broader context of Galactic structure studies and emphasize the 
absolute need of accurate spectroscopic stellar parameters in order 
to fully exploit the enormous amount and quality of data emerging 
from large-scale stellar surveys. 



2 OBSERVED STARS, MODEL ATMOSPHERES AND 
STELLAR PARAMETERS 



The observed dataset was taken from dRuchti et al 



comprises the core of the sample used in lRuchti et al 



2011) and it 
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after Paper 1) as well. These stars were originally selected for 
high-resolution observations based on d ata obtained by the R a- 
dial Velocity Experiment Survey (RAVE^ Isteinmetz et al l2006h in 
order to study the metal-poor thick disk of the Milky Way. High 
signal-to-noise spectra for these stars were obtained using high- 
resolution echelle spectrographs (R > 30000 and S/N ~ 100 
per pixel) at several facilities around the world. For the present 
work, the final sample contains 254 stars in the metallicity range 
-3.0 < [Fe/H] < -0.5. With respect to the full sample, we have 
excluded horizontal branch candidates because the grid of stellar 
models used here extends only up to the tip of the red giant branch 
(RGB). Also, stars with logg < 0.5 have not been considered to 
avoid large extrapolations fr om the grid of N LTE corrections that 
extends down to logg = 1 dLind et al]|2012if1 Stellar parameters 
were determined using two different techniques, as described in de- 
tail in Paper 1 and summarized below. 

First, we apply the classical method of ID LTE excitation- 
ionization equilibrium of Fe 1 and Fe n, which relies on the min- 
imization of the abundance trends with reduced line equivalent 
widths and lower level excitation potential. T e g and microturbu- 
lence are optimized to establish excitation balance, and logg and 
[Fe/H] are found by simultaneously requiring ionization balance. 
Notwithstanding the physical limitations due to the assumptions of 
ID LTE, the method is still in widespread use today because of its 
simplicity and easiness of implementation. 

In the second approach, substantial efforts were made to im- 
prove the accuracy of basic stellar parameters in the attempt to 
minimize the above-mentioned systematic errors. Firstly, effec- 
tive temperatures were determined from different methods avail- 
able for cool stars (such as the Balmer lines, semi-empirically cali- 
brated photometry). The final, optimal, T eff is that derived from the 
weighted averaging of different methods. In a second step, we re- 
computed gravities and metallicites relaxing the LTE assumption 



computed completely in NLTE. Such a code is currently being developed by 
one of the co-authors and will be applied to a large dataset of low-resolution 
spectra. 

2 The NL TE abundance cor rections for the Fe I lines are available for the 
full grid of lLind etafl i2012l) under |www . inspect- stars . net| 
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Figure 1. Comparison of stellar parameters derived using the LTE-Fe method and the NLTE-Opt stellar parameters for our sample of metal-poor stars. 
Differences in temperature, gravity, and metallicity show a large systematic increase with decreasing metallicity. See Paper 1 for full details. 



and adopting the new T cft scale. The NLTE corrections were ap- 
plied to Fe i abundances using the new wrapper code to MOOG. 
The typical errors are a(T c jf) as 60- 100 K, <x(logg) a 0.1 dex, and 
cr([Fe/H]) « 0.1 dex. Note that these are the absolute errors (sys- 
tematic + internal) on the accuracy of our stellar parameters. Sim- 
ilar errors bars are typically quoted in standard techniques. How- 
ever, usually those represent the internal precision of a method only. 
The results obtained by the LTE-Fe method are affected by much 
larger systematic errors. In the Paper 1, we contrasted our results 
for stars with surface gravities and temperatures determined by 
independent techniques, i.e., from accurate Hipparcos parallaxes 
and interferometry. We found that while the LTE-Fe surface grav- 
ities were systematically too low by -0.32 ± 0.39 dex, the NLTE- 
Opt surface gravities were in agreement with astrometric values to 
0.02 ± 0. 1 1 dex. This supports the accuracy of the new method. 

The spectroscopic stellar parameters obtained using both 
methods described above are compared in Fig. [TJ Hereafter, these 
two sets of stellar parameters are referred to as 'Plte.Fc (LTE-Fe 
method) and fNLTE.opt (NLTE-Opt method). The major conclusion 
to be drawn from these plots is that the classical method of ID LTE 
excitation-ionization balance of Fe leads to systematic errors in all 
three basic stellar parameters: 

• 7W (Fig.[TJ top panel): the LTE-Fe T s g values are too low by 



100 to 400 K for all stars below [Fe/H] 1; for dwarfs above this 

metallicity, the effective temperatures are usually over-estimated. 
Note that the LTE-Fe r e ff values for very cool giants appear to be 
reasonably accurate for [Fe/H] > -1.5. 

• logg (Fig. [TJ middle panel): surface gravities are usually 
under-estimated. The largest effects is seen for warm giants, T c ff > 

4500 K, which are biased by up to 1.5 dex at [Fe/H] 3 and by 0.5 

dex at [Fe/H] ~ -1.5; for metal-poor dwarfs, gravities are under- 
estimated by ~ 0.4 dex. For stars with [Fe/H] > -1, the errors are 
within 0.2 dex. 

Fe/H (Fig. [TJ bottom panel): the error is nearly a linear function 
of metallicity itself, such that for dwarfs and giants with [Fe/H] 

> -1 the error in [Fe/H] is 0. ldex. The error changes sign and 

increases to +0.2 dex at [Fe/H] 1.5, while at [Fe/H] 3 the 

LTE metallicity is too low by 0.5 dex for giants and 0.3 dex for 
dwarfs. 

Is it seen that the systematic effects are metallicity- and gravity- 
dependent, which is reflected in a notable spread of Tnlte.opi ~~ 
r LTE ,Fc and log g N LTE,opi - log gLTEFe at a given log g and r cff . 

Fig.[2]compares the location of the stars for the parameter sets 
^lte.Fc (left panel) and fWraopt (right panel) on the T cff - logg 
plane together with refere nce stellar evolutionary tracks computed 
using the GARSTEC code dWeiss & Sch lattl 2008, see next section 
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for a description of the grid of stellar models). Tracks are identical 
in the left and right plots. Symbol sizes are proportional to [Fe/H] 
as indicated in the right figure. Even though the tracks are shown 
only for visual guidance, it is clear that the ^NLTE.opt parameters 
trace the morphology of the turn-off much better. Results for gi- 
ants are also quite encouraging. For each star with T cft < 5300 K, 
we have searched the entire grid of stellar models the find the best 
fit to ? > LTE.Fe an d '? > NLTE,opt parameters as defined by the minimum 
X 2 . In the case of 'Pi,TE,Fe parameters, the average minimum^ 2 for 
these stars is 1.9, whereas if 'Pnlte.opi parameters are used it is 0.4. 
Clearly, 'Pnlte.opi parameters lead to an improved overall agreement 
with stellar tracks. 



for [Fe/H] < 0.0 and 0.05 dex for [Fe/H] > 0.o[J The set of tracks 
on which we base the results presented in the following sections in- 
clude microscopic diffusion and a 0.4 dex tt-element enhancement 
for [Fe/H] < -0.6. A maximum depletion of metals due to gravi- 
tational settling of about 0.4 dex is obtained for stars typical of the 
turn-off in metal poor ([Fe/H] ~ -2) globular clusters. No addi- 
tional extra mixing to those mentioned is included in the models. A 
detailed discussion on this is of no relevance to this work because 
we are interested in assessing the effect of input spectroscopic pa- 
rameters in the determination of stellar properties, and these re- 
sults do not depend on the reference stellar tracks used. Models are 
evolved from pre-main sequence up to the tip of the RGB or the 
age t = 30 Gyr, whichever occurs first. 



3 EVOLUTIONARY TRACKS 

Stellar models have been computed with GARSTEC and detailed in- 
formation about numerical aspects of the code and th e input physics 
available can be found in IWeiss & Schiattil J2008h : here we just 
provide a very brief summary of the most relevant input physics 
included in the stellar models. 

Regarding microphysics, we h ave used the FreeEOS0 equa- 
tion of state dCassisi et ail l2003h and the 14 N(p,y) ' 5 reac- 
tion rate recommended by LUNA dMarta et ail l2008h . Convec- 
tion is accounted for by using th e standard mixing length the- 
ory dKippenhahn & Weigert|[l990l) and the mixing length param- 
eter q-mlt is taken from the calibration of a solar model (see be- 
low). Convective overshooti ng is modeled as a diffusive process 
fol lowing the formulation by dFrevtag et al.ll9961 and see also Eq.3 
in IWeiss & Schiattil 1200 8h and the free parameter / is fixed to 
0.02. This value represents a moderate amount of overshooting 
on the main sequence comparable to the canonical 0.2 Hi0- How- 
ever, overshooting has to be limited in small convective core|fl In 
GARSTEC this is now achieved by implementing a geometric cutoff 
where the extension of the overshooting region decreases quadrat- 
ically with the ratio between the size of th e convective core a nd 
the pressure scale height at its boundary (see lMagic et al.ll2010l for 
details). 

In order to fix the mixing length parameter and the reference 
composition of our models we have cal ibrated a standard solar 
model using the|G revesse & Sauvall dl998l) solar composition, from 
which the present-day surface (Z/X) Q = 0.0229 and log epe = 7.50. 
The latter defines [Fe/H] = in our models. From the calibration 
we obtain q-mlt = 1.811 and the initial solar helium abundance 
Y m = 0.26896 and metallicity Z INI = 0.01876. Due to the effects 
of microscopic diffusion, the initial solar composition corresponds 
to [Fe/H] = +0.06. 

For the present work, we have computed a dense grid of stel- 
lar models in the mass range 0.6 < M/M Q < 1.4 in steps of 
0.01 Mjf|. The initial composition of the models is computed as- 
suming a cosmic helium-to-metal enrichment AY/AZ = 1.4 an- 
chored to the initial solar composition. The initial Fe abundance in 
the grid spans the range -5.0 < [Fe/H] < +0.5 in steps of 0.1 dex 



3 http://freeeos.sourceforge.net 

4 Hp = P/gp is the pressure scale height, where P and p are the pressure 
and gas density respectively. 

5 For example, a solar model computed with the canonical amount of over- 
shooting of 0.2 Hp still has, at the present solar age, a convective core. This 
is ruled out from helioseismic results. 

6 The stellar evolution models can be provided upon request to the main 
author 



4 BAYESIAN ANALYSIS 

Bayesian inference offers a powerful way to characterize a sys- 
tem by allowing the incorporation of prior knowledge into the sta- 
tistical analysis in an objective way. In the context of our prob- 
lem, i.e. determination of fundamental stellar parameters, Bayesian 
methods have been discussed and developed by var ious authors. 
Of par ticul ar pedagogical interest are the articles by IPont & Eyed 
d2004l) and[j0rgensen & Lindegrenl d2005l) and more recent and so- 
phisticated a nalysis with applications t o large stellar s urvey s can 
be found in iBurnett & Binnevl J2010h . iBurnett et ail d201ll) and 
ICasagrande et al. ( 2011 ). Our implementation of Bayesian infer- 
ence of basic stellar parameters is described in the next section. 



4.1 The Method 

In Bayesian analysis, the stellar parameters we want to determine, 
e.g mass and age, are random variables for which the probability 
distribution function (PDF) needs to be found. In particular, the 
posterior PDF, that is the probability of a model being correct given 
the observed data and the prior knowledge, is expressed as 



p(M,T, [Fe/H] | r eff , D 



■) 



p(M, t, [Fe/H]) 



x L(T el 



| M, t, [Fe/H]) 



(1) 



where the function p(M, t, [Fe/H]) on the right-hand side rep- 
resents the prior probability of the stellar parameters and 
L(r eff obs , ... | M, t, [Fe/H]) the likelihood, i.e. the probability of the 
observed parameters given a model characterized by a set of ob- 
servables D,(M, r, Fe/H]). Assuming independence and Gaussian 
errors in the observed data, the likelihood is 
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In our case, D consists of three parameters: r eff , logg, and [Fe/H]. 
Other parameters, like photometric information, can be readily im- 
plemented, but we do not consider them here. 

The prior probability density of model parameters can be writ- 
ten as 

p(T, M, [Fe/H]) = P (t) x p([Fe/H] | t) x p(M | [Fe/H], t) (3) 



7 The initial [Fe/H] in the evolutionary tracks is 0.06 dex higher than 
quoted on the label. We keep this nomenclature in order to associate 
[Fe/H]= evolutionary tracks to those with initial composition equal to 
initial solar composition. 
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Figure 2. The location of the stars for the parameter sets 7 3 LTE,Fe (left panel) and Pnlte.opi (right panel) on the HRD. For comparison, the reference stellar 
evolutionary tracks of 0.8 and 1 M Q and metallicities indicated in the legend are overplotted. Symbol sizes are proportional to the metallicity of the observed 
stars in both panels as indicated in the right plot. 



For the prior on mass, we assume a standard Salpeter initial mass 
function (IMF) independent of age and metallicity. For metallicity 
we assume a flat prior, i.e. that at all ages all metallicities are a 
priori equally probable. This simple choice is justified in this work 
because we are not interested in studying age-metallicity relations 
for stellar populations and because sufficiently reliable estimates of 
metallicity are available. Finally, for the age we also assume a flat 
prior, i.e. that star formation rate (SFR) is constant, but we set a cut 
at T max = 15 Gyr. Therefore 



p(M) ~ in ' 
p(r) ~ 1/t„ 



m 6 [0.6, 1.4] M 
t e [0, 15] Gyr. 



(4) 

(5) 



More complex priors can be elaborated, but our main goal here is to 
investigate the differences in masses, ages, and distances due to the 
input spectroscopic parameters rather than to characterize stellar 
populations. 

Finally, the proportionality constant in Eq. Q] 
is the normalization of the posterior PDF such that 
Jjfjf p(M, t, [Fe/H] | r eff , obs ,...)dMdTd[Fe/H] = 1, and it is 
of no relevance for the estimation of parameters and uncertainties. 

The PDF of a stellar parameter x that takes the values x tI along 
evolutionary tracks is given by 



p(x I r cff , obs , ...) = fffS(x- x a ) 

p(M, t, [Fe/H] | r 6fr , 0bs , ...)dMdTd[Fe/H] 

where the integration is over all evolutionary tracks. 



(6) 



4.2 General properties 

A key advantage of Bayesian inference is that, as expressed before, 
it allows using prior knowledge about the problem being consid- 
ered. Priors are particularly important when a quantity of interest, 
affected by the prior, has noticeable changes within the range of 
observational uncertainties where the likelihood doe s not change 
much. A typical example, ori ginally discussed by | Pont & Even 
and further expanded by |j0rgensen & Lindegrenl J2005h . is 
the different evolutionary speed of stellar models that populate the 
same region of observed parameters. When this is not accounted 
for, as in methods purely based on maximum likelihood, the results 



are biased towards faster evolutionary phases, a bias known as the 
'terminal-age bias' jPont & Everll2004l) . 

To illustrate here the importance of using information not 
contained in the likelihood, we present a similar example. 
Let us assume a star has the following observed parameters: 
T cff = 6850 ± 50 K and logg = 4.05 ± 0.05 dex. For simplicity, we 
exclude metallicity in this example by fixing it to [Fe/H] = -1.0. In 
the top panel of Fig.[3]we show three evolutionary tracks of models 
with different masses overlapping in the range of the observed pa- 
rameters. While the 1.1 M Q model is entering the subgiant branch 
(SGB), the 1.2 M Q is on the main sequence (MS) and the 1.3 Ms 
model is still on the pre-main sequence (PMS). 

First, we perform a standard maximum likelihood analysis 
(hereafter, ML) to determine the stellar mass and age based on 
Monte-Carlo sampling. We create synthetic distributions of T eff and 
log g assuming Gaussian distributions of errors and search for the 
best-fit model in all evolutionary tracks of the assumed metallicity. 
The resulting mass and age distributions are shown in the middle 
and bottom panels of Fig. [3] respectively. The mass distribution 
shows a bimodality with a secondary peak and a tail extending to- 
wards larger masses (M > 1.2 M Q ). The age distribution shows a 
dominant peak at very low ages (6 Myr < r < 10 Myr) clearly seen 
in the inset that zooms in in this age range. There are two other local 
maxima at about 4.3 Gyr and 4.8 Gyr, but with much lower relative 
probabilities than the solution at very low age. Thus the ML results 
indicate a strong preference for the PMS solution for the observed 
parameters. 

While it cannot be completely ruled out that the star is on 
the PMS, considerations of the shortness of this evolutionary phase 
render this possibility, a priori, not very likely. A 'brute force' way 
of taking this into consideration is to directly eliminate PMS mod- 
els from the analysi^. We impose this restriction by removing all 
stellar models with ages t < 20 Myr. In this way, while the best-fit 
model is still determined purely from its ML, we are already us- 
ing prior information by assuming the chance of observing a PMS 



8 This has been done, for example, in asteroseismic studies iGai et al] 
1201 ll) because PMS models have global oscillation properties that can be 
degenerate with those of MS counterparts. 
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Figure 3. Top panel: degeneracy of evolutionary tracks on the r c ff - logg 
plane at the location of our fiducial star (see text). Middle panel: relative 
probability of the mass distribution functions obtained using pure maximum 
likelihood (ML - blue dotted line), maximum likelihood with suppressed 
pre-main sequence models (ML no-PMS - red dashed line), and Bayesian 
analysis (Bayes - black solid line). Bottom panel: same as above but for 
the age distribution functions; the inset shows in detail the physically- 
disfavoured narrow peak corresponding to the pre-main sequence solution 
(see text for details). 



star is negligible (actually null, although this also is a wrong as- 
sumption!). Results for this exercise are labelled 'ML no-PMS' in 
the middle and bottom panels of Fig. [3] In the mass distribution, 
the tail at high masses vanished since its origin was purely from the 
PMS models. The same is true for the ages: the unrealistically large 
peak at ~ 6 Myr disappeared, and the peaks at 4.3 Gyr and 4.8 Gyr 
are now the dominant features, but their relative weight is still the 
same as in the pure ML analysis and the older peak, associated with 
stars entering the SGB still dominates. 

The results of the Bayesian analysis for our fiducial star are 



also shown in Fig. [3] (middle and bottom panels, solid black line). 
As in the case of the ML case, three peaks are present but the maxi- 
mum probability in each of them reflects information not contained 
in the ML analysis. The PMS solution, although possible, is now 
strongly disfavoured, whereas the MS solution is now the most 
likely one and it outweighs the SGB solution. The Bayesian mass 
distribution (middle panel) is similar to the ML no-PMS case. The 
vanishing contribution of the PMS is now easy to understand from 
the eq.[6] The probability for a given mass is obtained by marginal- 
izing the posterior PDF over all possible ages. The stellar tracks 
with high masses, which contribute in the observational box with 
the PMS, are strongly disfavoured because of the short duration of 
this phase. However, the PMS solution still has a non-negligible 
probability in the age distribution. The widths of the peaks reflect 
evolutionary speed in each phase, with PMS evolution being about 
100 times faster than late MS evolution and early SGB a few times 
faster than late MS. 

There is one important aspect related to the interpetation of 
the height of the peaks in the Bayesian age PDF. Strictly speak- 
ing, when determining the probability of a given age p(r) using 
Equation [6] there is no information on the speed of stellar evo- 
lution because we are lookin g at a time snapshot. As stated by 
|j0rgensen & Lindegrenl d2005l) . the relative heights of the peaks in 
the age distribution are determined by how closely packed evolu- 
tionary tracks of different masses are around the observational error 
box (i.e., the region over which likelihood is significant and varies 
little). Where variations in stellar mass introduce small changes on 
the location of a stellar model (for the same age), the integrand in 
Equation [6] will vary slowly when integrating in mass. Thus, the 
interpretation of the height of the peaks as depending on the evolu- 
tionary speed of stars is not straightforward. It only follows if slow 
evolutionary speeds are associated with regions on the observa- 
tional plane where small changes in mass produce small changes in 
the observational quantities the determine the likelihood. As stated 
above, the unique direct manifestation of a slowly evolving phase 
in the age PDF is the width of the peaks in the distribution. 



4.3 Examples 

Fig. [4] shows in black thick lines the PDFs of stellar parame- 
ters of two stars from our sample. On the left, Jl 12344.9-290327 
has [Fe/H] = -0.46 ± 0.10, T, B = 6340 ± 70 K, and logg = 
4.12 ± 0.10, i.e. it is a moderately metal-poor turn-off star. In 
this case, mass, age and absolute Ks magnitud^l have PDFs that 
lead to well determined c entral values and uncertainties. Following 
iBurnett & Binnevl d2010l) . we compute stellar parameters and their 
uncertainty directly from the different moments of the PDFs. In this 
case, we obtain t = 4.78 ± 0.88 Gyr, M = 1.108 ± 0.071 M Q , and 
M Ks = 2.371 ± 0.279. 

As a second case we have selected Jl 15337.3-020036, with 
[Fe/H] = -2.03 ±0.10, T eff = 5287 ±68 K, and logg = 2.98 ±0.10, 
a star on the lower part of the RGB. As it is well known, evolution- 
ary tracks of stars with different masses and similar metallicity are 
almost degenerate along the RGB in the T e g - logg plane (Fig. [2J. 
This degeneracy precludes the determination of its mass and age 
as seen from the PDFs in Fig. [4] The cut in the age distribution at 
15 Gyr simply reflects our age prior, and this corresponds to the 



9 Hereafter, we refer to 2MASS Ks magnitudes only, because they are used 
to determine distances. 
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Figure 4. Probability distribution functions (PDFs) for age, mass, and ab- 
solute magnitude in the Ks band. On the left results are shown for a turn-off 
star and, on the right for a star on the red giant branch. In the latter case, 
neither mass or age can be determined. The absolute magnitude (luminos- 
ity), on the other hand, can be derived precisely in both cases. Black solid 
lines: Salpeter IMF as a prior in the mass distribution; red dashed lines: flat 
IFM. 



mass cut at ~ 0.75 M Q . The shapes of the PDFs are almost com- 
pletely determined by the adopted priors, i.e., IMF and SFR, and 
by the fact that more massive stars evolve faster along the RGB 
than lower mass stars. No actual information about the true stellar 
mass and age can be obtained in this case. In fact, a pure likelihood 
analysis leads to age and mass distributions that are basically flat. 
On the other hand, the absolute magnitude of this star is still well- 
defined, and in fact yields M Ks = -0. 155±0.289, the same precision 
as for the turn-off star. This conclusion remains valid for all stars 
along the RGB, the absolute magnitude is a well-defined quantity 
even if the underlying evolutionary tracks are highly degenerate, 
and we can robustly estimate stellar distances independently of the 
evolutionary stage of the stars. 

Finally, we have tested the impact of the IMF prior by consid- 
ering also a flat IMF. Results for the two examples just discussed 
are shown in Figure |4] in red dashed line. Clearly, the impact of 
changing the IMF is minimal in those cases where stellar parame- 
ters are well defined. On the contrary, the large impact of priors in 
the age and mass distributions for the giant star are evidence that 
no actual information can be extracted in these cases: the posterior 
PDF is largely determined by the priors. 



5 RESULTS 

Using the procedure described in Sect. I4.1I and the two sets, 
fuE.Fe and '? > NLTE,opt of stellar parameters (Section |2}, we com- 
puted ages, masses and absolute magnitudes for the stars in our 
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Figure 5. Comparison of ages determined for the sample of nearby stars of 
iLachaume et alHl999l) . 
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Figure 6. Comparison of distance determination methods. We show dis- 
tances o btained with our Ba yesian scheme and those by the RAVE collab- 
oration iZwitter et al.ll2O10l) . Note that, for this test only, the input stellar 
parameters (r e (f, logg, [Fe/H]) are from the RAVE DR3 catalogue for both 
distance determinations. 



sample. To determine the distances, we used the observed 2MASS 
Ks magnitudes and the absolute Ks magnitudes obtained from 
the Bayesian analysis. The synthetic mag nitudes were computed 
from the MARCS grid of synthetic spectra dGustafsson et al.l2008h 
assuming plane-parallel radiative transfer for log g > 3 and 
spherically-symmetric mode otherwise. The synthetic Ks magni- 
tudes have been comput ed adopting the absolute calibrations and 
zero-points discussed in ICasagrande et all |2010|). Reddening for 



each star was computed using [Schlegel et al. ( 1998) maps follow 



ing the iterative procedure described in lRuchti et al 



fcOllI) . 



5.1 Tests 



In order to assess the robustness of our Bayesian scheme, we se- 
lected two sub-sample s of stars with basic stel l ar par ameters from 
Lachaume et al. l fl999T) and from lZwitter et al. For the for- 

mer, we determine age s and compare them with the results from 
Lachaume et al. I( fl99i) . For the latter, we focus on the distances. 
For stellar ages, the comparison only comprises stars for 
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which lEachaume et ail i 19991) derived ages from isochrone fittinj*^ 
in order to avoid biases derived from comparing ages obtained with 
methods not based on stellar evolutionary tracks. The very good 
agreement (Fig. [5} implies that when the observational uncertain- 
ties are small, 1% in T e g and 0.05 mag in M v , maximum likeli- 
hood methods return similar results to the Bayesian analysis and 
ag e determinations are robust. A similar agreement was reported 
by |j0rgensen & Linde gren (2005j). The uncertai nties in the ages 
we de rive are generally smaller compared both to lLachaume et al.l 
d 19991) and |j0rgensen & Lindegrenl ( l2005l) . 

We have derived stellar dista nces using the input stellar pa- 
rameters from lZwitter et aillfeOloH'I for the 58 stars common with 
our sa mple, and show the results in comparison with RAVE dis- 
tances dz witter et al.l20l6l) in Fig. [6] Results are an excellent agree- 
ment. Therefore, we conclude that, provided the same input stel- 
lar parameters are used, our Bayesian s cheme and the isochrone 
matching method bv lZwitter et al.l | |2010|) yield very similar results, 
with no obvious systemat ic offsets between methods. This is in line 
with lBurnett et al.l d201 ll) . who also applied a (more sophisticated) 
Bayesian analysis to the RAVE stellar sample and found that, al- 
though so me systematic differ ences were found, the overall agree- 
ment with lZwitter et ail ( 120100 was good. 



5.2 Ages and Masses 

Having verified the reliability of the Bayesian method, we proceed 
with the analysis of our stellar sample. The degeneracy of stellar 
tracks discussed in Sec. |4.3| limits the ability to derive masses and 
ages on the RGB. For this reason, we derived masses and ages for 
stars with logg > 3.6, that is, for 59 stars only. 

In Fig. |7] we compare masses and ages derived from the two 
sets 'PuE.fe and PNLTE.o P t °f input stellar parameters. Stellar masses 
derived from the PNLTE.opt parameters are systematically lower than 
those derived from Plte.Fo by about 0.05 - 0.1 M Q . The primary 
reason for the difference lies in the larger gravities obtained with 
the PNL.TE.opt method. Errors of the masses are comparable in both 
cases, and range between 0.03 M G up to 0.10 M Q . Nevertheless, the 
systematic shift to lower masses caused by the use of more accurate 
PNLTE.opt stellar parameters can not be neglected. 

For stellar ages, the behavior is opposite, as expected. The 
results based on the PNLTE.o P t parameters are higher by 10 to 30 
percent. However, not in all cases a smaller mass implies a larger 
age. Close inspection of both panels in Fig. [7] reveals only two 
stars with larger PNLTE,o P t masses, whereas six of them have smaller 
ages derived from PNLTE.opt parameters. Typical uncertainties of our 
Bayesian ages are of the order of 20% and, except for two cases, 
range between 10% and 30%. We conclude that the LTE-Fe method 
underestimates stellar ages by 1 - 4 Gyr, and the errors are more 
pronounced for hotter and more metal-poor stars. 



In addition, stars hotter than ~ 7500 K were excluded from the com- 
parison because the grid of stellar evolution models is not complete above 
those temperatures and because we use the MARCS model atmospheres 
iGustafsson et alj2008l) . which are available for 7/ cff < 8000 K only. 
" We used the RAVE distances obtained for stars with repeated obser- 
vations and Padova isochrones. This combination should provide relia ble 
estimations according to the note added in proof in lZwitteret"ai]feO10l) . 
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Figure 7. Masses and ages for our sample of stars derived with the NLTE- 
Opt and LTE-Fe basic stellar parameters. Only stars with log g > 3.6 are 
included (see text). 

5.3 Distances 

The main goal of this section is to show the dramatic impact that 
using the physically realistic PNLTE.opt stellar parameters has on the 
determination of stellar distances. 

As summarized in Sect.[2] there are large differences between 
the 'Plte.Fc and PNLTE.opt stellar parameters sets. The top panel of 
Fig. H] shows the change in the absolute magnitude Mk s as a func- 
tion of the log gNLTE.opt _ log gLTE,Fe differences. The strong linear 
correlation simply reflects the fact that changes in log g are accom- 
panied by a reciprocal change in luminosity due to their inverse 
functional dependence on stellar radius. Clearly, the PNLTE.opt stel- 
lar parameters indicate that all our stars, with few exceptions, are 
fainter than what would be inferred from using the LTE-Fe param- 
eters. Furthermore, if, for the sake of discussion, we neglect the ef- 
fect of reddening, then distance scales exactly as the square root of 
the luminosity of the star. Consequently, distances derived from the 
P LTE Fe set of parameters are overestimated, as illustrated in Fig. [8] 
(bottom panel) where we show the ratio of distances derived using 
pLTE,Fe and Pnlte,o P i stellar parameters. 

Similar information is provided in Fig. |9] in which the ratio 
of distances, corrected for reddening, is plotted as a function of 
loggLTE.Fc- This illustrates clearly how systematic errors in distance 
determinations impact different stellar types. Stellar metallicities 
are encoded in the symbol sizes and colours, as indicated in the leg- 
end. Distances to metal-rich stars, [Fe/H] > -1, are, to first order, 
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Figure 8. Changes in the absolute Ks magnitudes and distances caused 
by using the improved set of input basic stellar parameters, T^NLTE.opi as a 
function of the change in log g, for our full sample. 
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Figure 9. Ratio of distances derived with PuE.Fe and !PNLTE,opt stellar pa- 
rameters as a function of surface gravity for our full sample. Symbol sizes 
and colours indicate the stellar metallicity. 



unbiased. On the other hand, for metal-poor ([Fe/H] < -1) dwarfs 
and sub-giants (log g > 3.6), LTE-Fe parameters imply a systematic 
distance error of about 10 - 40%. The problems are exacerbated at 
low metallicity and low gravity. Metal-poor giants suffer from the 
largest fractional distance biases of up to 70%. 

Fig. [10] shows the final distribution of distances determined 
using the fNLTE.opt and !Plte,Fc stellar parameters for our full sam- 
ple. This plot summarizes what has already been discussed above, 
and allows us to put the results into the more general perspec- 
tive of large-scale stellar surveys. In a magnitude-limited survey, 
for example, where more metal-rich unevolved stars dominate the 
nearby sample and metal-poor luminous giants are predominantly 
observed at large distances, classical full LTE-Fe analysi^ will 



It is important to remind the reader that LTE-Fe refers to the classical 



systematically overestimate distances, placing stars progressively 
further than they are. This has important consequences in all cases 
when stellar metallicities, positions, and kine matics are used for 
studies of Galactic structure and evol ution (Siebert et al. 
ICarollo et al l |2010| : Iwilson et al.l 1201 ll ; iBensbv et all 1201 it) . We 
touch upon this in Sect. [6] 



5.4 Comparison with other studies: RAVE stellar parameters 

Since our sample of stars was drawn from the internal RAVE cat- 
alogue, it is very important to investigate how our results, i.e. r e ff, 
logg, [Fe/H], and distances, compare with those obtained by the 
RAVE collab oration. In what foll ows, we will compare our data 
with those by IZ witter etal] ( |2010| , R with which we have 58 stars 
in common. We highlight the key methodological differences be- 
low. 

• The quality of observations. Our parameters rely on high- 
resolution data, R ~ 30000, with a very broad wavelength cover- 
age, 3500 - 9500 A; the RAVE spectra have R ~ 7500 and cover 
only a limited wavelength range between 8410 and 8795 A. 

• The physics of stellar atmospheres and radiative transfer mod- 
els. In particular, our analysis relies on NLTE radiative transfer, 
while the RAVE spectroscopic pipeline operates fully on ID LTE 
synthetic spectra. 

• The different approaches to determine masses and ages, and 
stellar e volution mode l s. GAR STEC tracks were used in this work, 
whereas IZ witter et al. I ( 1201 Ol) adopted Padova isochrones. Also, 
the approach followed in that study is esentially a nearest neigh- 
bor technique, however with the important difference that the 
isochrones were resampled to a uniform spacing in age to account 
for the variable speed of stellar evolution, and a weighing scheme 
to include, for example, the impact of the initial mass function. In 
this respect, their analysi s also makes use of priors. 

• IZ witter et al] fcOlOl) neglected reddening in their distance de- 
termin ations. We used the reddening values following iRuchti et al .1 
J20l ]]); also comparison with the interstellar NaD lines in our high- 
resolution spectra indicates that reddening is not negligible for 
many stars. 

• The uncertaintie s on our input stellar parameters are typically 
smaller than those of IZ witter et al.l d2010h . Their values are ~ 0.2 
dex in [Fe/H] and logg and range between ~ 100 to 200 K in most 
cases for r cff , whereas our standard estimates are 60- 100 K in r eff , 
and 0.1 dex in [Fe/H] and logg. 

The differences between ? 3 NLTE,opt parameters and those from 
the RAVE DR3 pipeline are shown in Fig. QT] Size and colour of 
the symbols indicate metallicity and gravity (we used RAVE DR3 
values for reference). We also show the differences fNLTE.opt and 
^LTE.Fe for the same stars in grey empty squares. RAVE DR3 pa- 
rameters show a very large spread with respect to ^NLTE.opt • In fact, 
the spread is much larger than with respect to the 'PrjE.Fe param- 
eters. Clearly it is not possible to single out specific factors which 



spectroscopic method of LTE excitation-ionization balance of Fe, which is 
very efficient because T c fj, logg, and [Fe/H] can be all determined simulta- 
neously; 

13 ISiebert et al] feplll) have presented the revised RAVE DR3 pipeline. 
However r c u and l og g for the 58 stars we use in this comparison are ex- 
actly the same as in IZwitter et allfeOld) ; only [Fe/H] differ by about 0.1 to 
0.2 dex, and this has a very minor im pact on the de t ermin ation of distances. 
We thus prefer to use all data from 
used by RAVE to obtain distances. 



Zwitter et al.l fcOlOl) . which has been 
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Figure 10. Comparison of heliocentric distances derived for all stars in our 
sample using !pLTE,Fe and 7 3 NLTE,opi sets of basic stellar parameters. In virtu- 
ally all cases, ?\,TE,Fe parameters lead to large overestimation of distances. 
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Figure 11. Differences in basic stellar parameters for the stars in com- 
mon between our sample and the RAVE distances catalogue (Zwitter et al. 
l201tj) . Colours and symbol types are assigned according to logg values 
from the DR3 pipeline. Symbol sizes are indicative of the DR3 metallic- 
ity, as indicated in the legend. The differences between (NLTE,Opt) and 
(LTE.Fe) parameters for the same stars are shown in grey empty squares. 




cause the spread; however, some conclusions can be drawn. In most 
cases, the RAVE DR3 logg values are clearly underestimated for 
giants with A log g ( NLTE,o P t-RAVE DR3) up to ~ 1 dex and one case 
as large as 1.3 dex. This effect is more pronounced for the metal- 
poor stars, in qualitative agreement with the differences between 
the 'Pnlte.opi-'^'lte.Fc parameters for the same stars. The same is true 
for subgiants; in most cases the DR3 gravities are underestimated. 
For dwarfs, on the other hand, we find DR3 gravitie s are overes- 
timated by as much as 0.5 dex. iBurnett et alj d20lTh arrive at the 
same conclusion comparing DR3 gravities with stellar models. 

As a consequence of the differences in stellar parameters, par- 
ticularly changes in log g, the distances we obtain from the 'Pnlte.opi 
parameters are systematically smaller than the RAVE DR3 data 
jZ witter et al.l20irilsiebert et al.l201 lb . For the most distant metal- 
poor giants, the offset is of the order 50 percent. This is illustrated in 
Figure [L2l which shows that beyond 1.5 kpc all ^NLTE.opt dist ances 
but one are smaller than those derived in lZwitter et alj (HoiOh . For 
a few dwarfs in the immediate solar neighborhood, the results are 
opposite, as expected from the sign of the logg corrections. The 
overall trend is very clear and the qualitative similarity with the 
difference between 'Pnlte.opi and ^LTE.Fe distances (Fig.llOt is strik- 
ing. 

These results provide alarming evidence that the RAVE DR3 
stellar parameters are affected by the systematic error caused by 
the ID LTE assumption in spectroscopic parameter determina- 
tions. Test calculations for t he Call near-IR triplet us ing the mean 
3D model atmospheres from lBergemann et alj d2012l) indicate that 
these spectral lines are strongly affected by the ID - (3D) differ- 
ences, while the near-IR Fe I lines are mostly affected by NLTE. 
Fig.EHillustrates these effects for the two well-studied metal-poor 
stars, a subgiant HD 140283 (top plot) and a warm turn-off star HD 
84937 (bottom plot) with [Fe/H] = -2.4 and -2 de x, respectively. 
The spectra are from [Bergeman n & Gehrenl d2008b and shown in 
the top panels in their original resolution, R ~ 40 000. In the bot- 
tom panels the spectra have been degraded to the quality of the 
RAVE spectra, i.e. R ~ 7 500 and a typical S/N ~ 50 Jzwitter et ail 
l2008l Fig. 8). The (3D) NL TE synthetic spectra were computed as 



Bergemann et alj d2012h , but applying LTE radiative transfer to 



Calo The ID LTE radiative transfer models under-estimate the 
strength of the Call lines (the effect of (3D) structure) and over- 
estimate that of the Fe I lines (the effect of NLTE). Also, the near- 
IR lines of a-elements (Si I, Mg I, a nd Ti I), which alon g with Fe 
I lines do minate the RAVE spe ctra (Zwitt er et al. |l2008l) all form 
in NL TE dZhao & Gehrenll200rl IShi et al.|[2oTTl : iBergemann etai] 
I2OI2I) and it is unlikely that the NLTE effects cancel out to produce 
a spectrum close to LTE. The effects are clearly present irrespec- 
tively of the spectral resolution. 

We close this secti on with a few comm ents related to the two 
different validations that lz witter et al. l d2oToh presented for distance 
determination from RAVE data. The first one is a comparison with 
Hipparcos stars present in their catalogue with trigonometric errors 
smaller than 20% and spectroscopic distanc e errors small e r than 
32%. These results are shown in Figure 3 in IZ witter et al. 1 (12011 
and a large spread is present in the comparison. Even if no system- 
atic difference was found, as noted by the authors, almost all these 
stars dwarfs with a solar-like metallicity ([M/H]=-0.12 ± 0.25), 
exactly the case where we do not expect large NLTE effects. The 
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Figure 12. Comparison of distances derived using pNLTE.opt stellar parame- 
ters and those from the RAVE DR3 catalogue for stars common in both sam- 
ples. The trend strongly resembles that of the comparison between pLTE.Fe 
and fNLTE.opt distances (Fig.[T0l 



14 iMashonkina et alj 120071) found that NLTE effects in the near-IR Ca lines 
(3d-3p multiplet) are also important, and they may further strengthen the 
lines at low metallicity. 



Ages, masses, and distances of cool stars 1 1 



HD 140283 




8550 8600 
wavelength [A] 

HD 84937 



8650 



1.0 
0.8 
0.6 
0.4 
1.0 
0.9 
0.8 
0.7 
0.6 




Ca 



Cal 



-I 



I am- 40000; 

- LTE 1 D 

- NLTE,<3D> 




8450 8500 



8550 8600 
wavelength [A] 



8650 



8700 



Figure 13. Synthetic ID LTE (blue) and <3D> NLTE (red) spectra of the 
metal-poor stars HD 140283 (top plot) and HD 84937 (bottom plot). The 
observed spectra (black points) are from Bergemann & Gehren (2008) in 
their original resolution, R ~ 40 000 (top panel in each plot), and degraded 
to the resolution of the RAVE spectra R ~ 7 500 (bottom panel in each plot). 
The ID LTE radiative transfer models under-estimate the strength of the Ca 
II lines and over-estimate that of the Fe I lines (shown in the inset). See text 
for more details. 



second test consisted in comparing about 14 stars in 7 open clus- 
ters with k nown distances. A gain, all these clusters have solar-like 
metallicity jDias et al.ll2002l) and we do not expect NLTE effects to 
play a centra l role. In any case, we note that distances derived by 
IZwitter et al.l d2O10L see their Table 2) for the cluster stars deviate 
from the cluster distances by anything be t ween -20 to +40%. These 
stars were reanalyzed by [Burnett et al.l (1201 ll see their Table 3) 
and results show a similar spread, although in a few cases distance 
determ inations are also quite different from those in Iz witter et al.l 
J2010h . It is simply not possible to try to establish confidence in 
spectroscopically derived parameters, particularly at low metallici- 
ties and gravities, by extrapolating confidence gained by means of 
tests performed in very different regimes, e.g. at solar metallicity. 



6 SUMMARY AND CONCLUSIONS 

This paper is the second in the series where we explore the influ- 
ence of different methods on the determination of basic and fun- 
damental stellar parameters for FGKM stars. The main goal is to 
understand to what extent our knowledge about stars, their surface 
and interior structure, is biased by inferences derived from classi- 
cal model atmospheres with their simplifying assumptions (ID hy- 
drostatic and local thermodynamic equilibrium, LTE), which have 
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Figure 14. Heliocentric spatial distribution of stars in our sample. Top and 
bottom panels respectively show results obtained based on PlteSv an d 
^NLTE.opt parameters. Colours depict the stellar surface gravity as indicated 
in the plots. 



been used as a basis for spectroscopic characterization of stars over 
the past 70 years. 

In Paper 1, we developed a novel NLTE-Opt technique to de- 
termine accurate effective temperature, surface gravity, and metal- 
licity from a stellar spectrum and then applied the method to a 
large sample of stars selected from the RAVE survey. The prin- 
cipal components of the method are (i) the use of a robust effective 
temperature scale from the analysis of Balmer line wings and, (ii) 
the determination of surface gravities and metallicities by solving 
physically-realistic NLTE radiative transfer in Fe lines. This avoids 
substantial systematic uncertainties in parameters introduced by the 
commonly-used approximation of LTE. A comparison of the new 
NLTE-Opt stellar parameters to that obtained from the standard 
method of ID LTE excitation-ionization balance of Fe, LTE-Fe, 
revealed important differences ranging from -100 K to +400 (sub- 
giants, RGB) in r c ff, 0.1 to 1.5 dex in logg, and 0.05 to 0.5 dex in 
[Fe/H]. 

In this paper, we focus on the determination of evolution- 
ary parameters of stars, masses and ages, and their distances. We 
develop a robust Bayesian technique using state-of-the-art evolu- 
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tionar y tracks computed with the GARSTEC code dWeiss & Schiattil 
l2008h . The importance of the Bayesian sch eme in this context 
has been raised in several earlier studies dPont & Ever! |2004 
|j0rgensen & Lindegrenll2005l) . This method has a principal advan- 
tage over standard maximum likelihood (or minimum x 2 ) meth- 
ods in that it allows a simple way of introducing additional (prior) 
knowledge in the statistical analysis. Test comparisons of results 
obtained with our Bayesian scheme with other similar studies con- 
firm its reliability and accuracy. 

We then apply the Bayesian method to both sets of stellar pa- 
rameters, LTE-Fe and NLTE-Opt, for the sample of stars from our 
Paper 1. The stars were selected from the RAVE catalogue. They 
cover the full parameter space in Tea (3500 to 6500 K), logg (3500 
to 6500), and [Fe/H] (-3.0 to -0.5) and sample the local volume 
within ~ 10 kpc. Spectroscopic distances were computed for all 
these stars using the 2MASS photometry. Masses and ages were 
determined for the 59 dwarf and subgiant stars (log g > 3.6) only, 
the limit set by the degeneracy of stellar tracks of different mass 
along the RGB in the r cff - log g plane. 

From these calculations, the following conclusions can be 
drawn. 

• The classical method of LTE excitation-ionization balance of 
Fe, as a rule, underestimates Teff , log g, and [Fe/H] of FGKM stars. 
As demonstrated in Paper 1, for the most metal-poor giants, [Fe/H] 

3, the errors are of the order +0.5 dex in metallicity. This has 

implications for the metallicity distribution function of stars in dif- 
ferent Galactic components. For our stellar sample, which is repre- 
sentative of the thick-disk, the metal-poor end of the MDF steep- 
ens and is s kewed towards the more metal-rich domain (by A [Fe/H] 
— 1-0.3 dex; lBergemann et alj|2013l) . 

• LTE-Fe stellar parameters result in overestimated distances 
compared to the NLTE-Opt results, and the bias depends on the 
metallicity and evolutionary stage of a star. The changes in surface 
gravity are most critical for metal-poor RGB stars, for which the 
NLTE-Opt distances are smaller by up to 70%. The distances to 
metal-poor subgiant and tumoff stars are in error by a factor of 2. 
Only more metal-rich turn-off stars are barely affected by this bias. 
Fig.Elshows the spatial heliocentric distribution of the stars in our 
sample determined with Pue.Fc (top panel) and ^ViLTE-opt (bottom 
panel) parameters and makes evident the large impact of using inac- 
curate stellar parameters for deriving spectrophotometric distances, 
in particular for giants with log g < 2.0. The overall distribution 
of distances in our sample steepens towards smaller distances as 
shown in Fig. [15] As a consequence, significant errors may enter 
the determination of space velocities and orbital eccentricities, and 
differences in the population membership of stars (i.e. thick disk or 
halo) are expected. 

• C omparing our NLTE -Opt results with the RAVE DR3 dis- 
tances dz witter et ai]|2010h for the stars common in both samples 
we find a pattern similar to the comparison between NLTE-Opt and 
LTE-Fe results. This is not surprising: metal lines that dominate the 
near-IR RAVE spectra are known to be strongly affected by NLTE 
effects (Sect. |5.4b and our test calculations with mean 3D models 
also show that the lines are sensitive to the ID - (3D) atmospheric 
structure differences. As a result, more distant stars (predominantly 
metal-poor giants) have larger DR3 distances compared to our val- 
ues (Fig.[l2]l, by 10 - 50%. 

In this respect, it is very interesting to interpre t our results in 
the light of the findings by ISchonrich et al.l d201 ll) where the au- 
thors concluded, based on simple analytical considerations that dis- 
tances for SDSS/SEGUE metal-poor stars are systematically over- 
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Figure 15. Distribution of distances for our full sample showing the global 
impact of using the physically realistic !PNLTE,opt stellar parameters to esti- 
mate distances. 

estimated. Consistent with the arguments presented in that study, 
our calculations predict that the systematic error in log g and, thus, 
distance increases with decreasing metallicity; [Fe/H] is affected 
too, although at a lesser degree. Thus, there will be two effects. 
First, more metal-poor stars will appear more luminous and thus 
more distant, which may be mis-interpreted for a very metal-poor 
outer stellar population. This goes in the direction of explaining 
why this component would appear to be more metal poor than 
the inner halo component. Second, the systematic error in distance 
may propagate into the error in transverse velocity, causing the ef- 
fect of a star on a retr ograde orbit. This strengthens the claim by 
ISchonrich et al. I J201 1 ) that the available obs ervational evidence for 
a counter-rotating Galactic halo component dCarollo et al.ll2010l) is 
weak. 

• Masses and ages are likewise affected in a systematic way by 
the differences between NLTE-Opt and LTE-Fe, although there is 
not a one-to-one mapping of the change in mass to the change in 
age. In the LTE-Fe method, the majority of stars appear younger 
and slightly more massive. The difference with NLTE-Opt is of 
the order 0.05 - 0.1 M slm . The ages are inaccurate by anything 
between 10 and 30 percent. In this respect, it will be interesting 
to re-investigate the age-metallicity relation in the solar neighbor- 
hood, which is usually drawn from stellar samples with spectro- 
scopic metallicities and ages obtained by simple isochrone fitting. 

We emphasize that, contrary to a wide-spread impression, de- 
partures from LTE and ID hydrostatic equilibrium will affect low- 
and high-resolution spectra. The reason is that these effects change 
equivalent widths (or, equivalently, line indices), and consequently 
abundances inferred from a diagnostic feature. This can not be 
compensated by decreasing resolution of a spectrum, which merely 
alters the line profile shape without affecting the total energy ab- 
sorbed in a line. 

In summary, we have shown that the accuracy of the stellar pa- 
rameters (r eff , log g and [Fe/H]) is crucial to obtain reliable masses, 
ages, and spectroscopic distances for late-type stars. Our findings 
lead to the inevitable conclusion that using the adequate r cff scale 
and accounting for the deviations from LTE and ID hydrostatic 
equilibrium in stellar atmospheres is of fundamental importance if 
we want to achieve an accurate characterization of stellar popula- 
tions. Until now, most of the work in the community has gone into 
developing sophisticated statistical methods to extract the largest 
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possible amount of information from large-scale stellar surveys. 
Although these efforts are undoubtedly valuable, they will not cure 
the major problem affecting the final distributions of masses, ages, 
metallicities, distances, which arises because systematic errors in 
basic stellar parameters dominate. We thus advise particular cau- 
tion when applying stellar parameters obtained from ID LTE mod- 
els in studies of stellar populations. 
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